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Abstract 



' By making use of the potentials of the heat conduction equation the integral equations are derived 

. which determine the heat kernel for the Laplace operator — a^A in the case of compound media. 

> 

0^ ' In each of the media the parameter acquires a certain constant value. At the interface of the 

oo ■ 

' media the conditions are imposed which demand the continuity of the 'temperature' and the 'heat 



flows'. The integration in the equations is spread out only over the interface of the media. As a 

result the dimension of the initial problem is reduced by 1. The perturbation series for the integral 

equations derived are nothing else as the multiple scattering expansions for the relevant heat 
Oh! 

^ ■ kernels. Thus a rigorous derivation of these expansions is given. In the one dimensional case the 

^ ■ 

integral equations at hand are solved explicitly (Abel equations) and the exact expressions for the 
^ ' regarding heat kernels are obtained for diverse matching conditions. Derivation of the asymptotic 

expansion of the integrated heat kernel for a compound media is considered by making use of the 
perturbation series for the integral equations obtained. The method proposed is also applicable to 
the configurations when the same medium is divided, by a smooth compact surface, into internal 
and external regions, or when only the region inside (or outside) this surface is considered with 
appropriate boundary conditions. 
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I. INTRODUCTION 



The heat kernel techniquo^ i^i^i^i^ is widely used for constructing the quantum field theory 
in gravitational background and with allowance for nontrivial boundary conditions. Of 
a particular interest is the asymptotic expansion of the heat kernel it terms of evolution 
parameter for its small values. The coefficients of this expansion pertain to divergences 
and anomalies in the relevant quantum field theory models. Proceeding from this one can 
develop the renormalization procedure needed. 

For well posed spectral problems the heat kernel coefficients are expressed, in an polyno- 
mial way, through the local geometric characteristics of the manifold D and its boundary S. 
Not only the contributions of D and S are independent but the contributions of individual 
regions of D and 5* are also additive. 

The spectral problem is well posed for the goal of constructing the heat kernel if the 
second oder elliptic differential operator in question is close to the Laplace operator defined 
on a smooth manifold with a smooth boundary, if any. 

There are no universal methods for constructing the heat kernel and its asymptotic ex- 
pansion. The development of different approaches to this problem is the subject of many 
works (see, for example, the reviewsi^^jMi^ and references therein) . 

The initial definition of the heat kernel is the Green function of the heat-conduction 
equation with an elliptic operator under study. In many physical problems it is worth going 
from the differential equation, defining the solution to be found or the relevant Green's 
function, to the equivalent integral equation. In the dynamical evolution problems the 
integral equations manifestly show the reason-consequence relations governing the physical 
process under study. Reducing the problem to the integral equation, as a rule, allows one 
to develop the method of successive approximations (perturbation theory). 

Transforming the initial differential equation into the integral form is in the general case 
a nontrivial problem. A special attention should be paid here to incorporating the boundary 
conditions into the integral equation. When constructing the integral equations governing 
the Green's function of the heat equation we shall use the surface potentials for this equation. 
The volume potential and the potentials of single and double layers naturally arise in the 
theory of the Laplace equation. In this case they are referred to as the Newtonian (or 
electrostatic) potentials. This idea proved to be fruitful also in studies of the Helmholtz 
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equation describing, for example, steady harmonic oscillations, the wave equation, and the 
heat-conduction equation. The potentials are particular solutions of these homogeneous 
equations, and they are constructed in a universal way in terms of the fundamental (or 
elementary) solution of the initial equation. The potential technique has turned out to be 
effective both for consideration of the general properties of the equations under study (for 
example, the prove of the solution uniqueness) and for deriving particular solutions with 
given properties and for obtaining the Green's functions. 

This paper seeks to demonstrate the efficiency of using the heat potentials when con- 
structing the integral equations for the heat kernel (Green functions) at first in the case of 
manifolds with boundary. It is worth noting that integration in these equations is spread out 
over the boundary only. As a result the dimensionality of the initial problem is reduced by 
co-dimension of the boundary 5*. Further this approach is extended to the compound media, 
where the principal part of the differential operator has a discontinuity at the interface be- 
tween different media. A typical example here is the electrodynamics of continuous media.— 
The velocity of light has, in the general case, a jump discontinuity on the border between 
two media with different characteristics (for example, on the border between dielectric and 
vacuum). In the both media the Maxwell equations are well defined, and at the interface 
the matching conditions (or boundary conditions) should be satisfied. The concrete form of 
these conditions is determined by the physical content of the problem in question. In the 
same way the conduction of heat in compound media is treated.— As far as we are aware, 
the heat kernel for compound media is not investigated yet.— 

The layout of the paper is as follows. In Sec. |n] the essentials of the potential theory 
are recalled first for the Laplace equation (Newtonian potentials) and then for the heat- 
conduction equation (heat potentials). By making use of the heat potentials the integral 
equation for the Green function (heat kernel) is derived for a compact region of Euclidean 
space bounded by a smooth surface. The perturbative series for this equation is developed 
which is nothing else as the multiple scattering expansion for the heat kernel. Thus a rigorous 
derivation of this expansion is presented. The convenience to use here the Laplace transform 
is shown. In Sec. IIIII the integral equations are derived that determine the heat kernel for 
compound media. The efficiency of the approach proposed is demonstrated by deriving the 
t-small asymptotics of the first three terms of the perturbation series for the heat kernel in 
the case of compound media f Section II V|) . In Sec. |3the heat kernel on an infinite line is 
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constructed in an exact form for diverse matching conditions. In the Conclusion f Section I VI|) 
the obtained results are briefly summarized and the possibility of extending the approach 
proposed is discussed. In Appendix |X] the general conditions at the interface are found 
which result in self-adjoint boundary value problem for the Laplace operator considered in 
compound media. 



II. HEAT POTENTIALS 



In order to recall the basic facts from the potential theory, we first address the Laplace 
equation 

Au = (2.1) 

considered in the d-dimensional Euclidean space R^, which is divided by a smooth closed 
surface S into a compact internal domain and external one D^. On the surface S the 
relevant boundary conditions should be imposed, that depend on the physical content of the 
problem in question. For example, it may be a mathematical formulation of the electrostatic 
problem.— The surface S is supposed to posses the properties of smoothness needed. In the 
potential theoryS'SiiSiii this implies that S is the Lyapunov type surface. The points of 
Euclidean space are denoted by x,y, z, . . ., and r^y is the Euclidean distance between x 
and y. At any point x on there exists a unit normal Ux or n{x). For definiteness we chose 
the inward directed normal. It is possible because we are dealing with closed surfaces S. 

For a linear homogeneous differential equation the fundamental (or elementary) solution 
is defined, which is the Green function of this equation in an unbounded space. By definition, 
the fundamental solution obeys the initial equation with a 5-like source on its right-hand 
side. In the case of the Laplace equation in M.^ the fundamental solution is 

with the properties 

A, Ea{x; x') = A,,, Ea{x; x') = 6'^''\x, x'), Ed{x; x') = Ed{x'; x) . (2.3) 

The potentials for the Laplace equation are constructed by making use of the fundamental 
solution, namely, the volume potential 

m.) = -^[^dy, (2.4) 
4 



the single-layer potential 



V{x) = -^ f'^dSy, (2.5) 



Arc Js r^y 



and the potential of a double layer 

V,(.).-_L/^,„)^J_,.„. (2.6) 

These formulas are written for d = 3, and the following notations are used: dy and dSy are, 
respectively, the elements of the volume and of the surface at the point y, w{y), fi{y), and 
ui^y) are the densities of these potentials. It is convenient to consider the densities w{y), 
fi{y), and z/(?/) to be continuous functions. 

All three potentials are solutions of the Laplace equation ()2.1|) . namely, the volume poten- 
tial ()2.4|) is harmonic outside D, V{x) and W{x) are harmonic outside S .The single-layer 
potential V{x) is continuous everywhere in W^, specifically, on passing through S. The 
potential of a double layer W{x) has a discontinuity on S, namely: 

W,{x) = W{x)~^fi{x), 

We{x) = W{x) + ^ /i(x) , xeS. (2.7) 

Here W{x) is the value of the integral ()2.6p . when the point x belongs to S {W{x) is a 
continuous function for x varying along S), Wi{x) is the value of the double- layer potential 
()2.6|) . when the point x tends to S from D+, and We{x) is the same when the point x 
approaches at S from 

In what follows we shall frequently use the derivative along the normal to the surface S 
at the point y, which belongs to S. This derivative acts on the function the argument of 
which is the distance r^y between the points x and y, the point x being not obliged to lay 
on S. Simple calculation gives 

d , , 

ra;y = COSV3, (2.8) 



dUy 



where is the angle between the vector Vxyi which starts at x and ends at and the normal 
My. In the same way we have 

d . w/ N Q ( ^\ cos 09 , , 

/(r.,) = fir^y) cos^, = . (2.9) 



duy ' dny \r^yj rly 

Equations ()2.8|) and ()2.9|) are obviously valid for both inward and outward directed nor- 
mals. 
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The application of the potentials ()2.4|) - ()2.7|) for transforming the boundary problems 
for the Laplace equation ()2.H) to the integral equations can be found in many textbooks on 
mathematical physics.— 

Let us proceed to consideration of the heat-conduction equation 

^-a^Au = 0. (2.10) 
The fundamental (or elementary) solution to this equation in is 

Ed{x, t- x\ t') = 9{t - t') Ko{x, t- x\ t') , (2.11) 

where 



Ko{x, t; x', t') = ^ exp 

2a^y^r{t - t') 



xx' 



4a2(t-f) 



(2.12) 



The function KQ[x,t] x' ,t') (the propagator— or the heat kernel) obeys the homogeneous 
heat equations 

Q \ 

— - a^A.^\ Ko{x,t;x',t') = 0, 

(^-^ + a^A,^ Ko{x, t- x\ t') = , (2.13) 

and inhomogeneous initial condition 

KQ[x,t;x\-i!)^ b^'^\x-x), when t ^ i! . (2.14) 

This condition enables one to construct the solution of the Cauchy problem for the nonho- 
mogeneous heat equation 

^-a'A,) n(x,t) = /(x,t) (2.15) 
considered in an unbounded space 

u{x,t) = / dx' Ko{x,t; x' ,t') Uo{x') + d9 dx' KQ{x,t;x',9) f{x',9) , (2.16) 



t' 



where uo{x) = u{x,t = t'). In the classical mathematical physical the representation ()2.16p 
of the solution to the heat conduction equation ()2.15|) is known as the Poisson formula. The 
last term in (j2.16p can be considered as an analog of Eq. (|2.4|) defining the volume potential 
for the Laplace equation. 
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In a complete agreement with the definitions ()2.5p and ()2.6|) the heat surface potentials 
are introduced, namely, the simple-layer potential 



2 



V{x, t)=a' de dSy Ko{x, t; y, 9) iy{y, 6) (2.17) 
Jo Js 

and the potential of a double layer 

W{x,t) = de dSy—^{x,t;y,e)fi{y,e) 

•-'0 ^ S y 

^^j—Q^ J dSy r cos (f Ko{x,t;y,e fi{y, 6) , (2.18) 

where the functions uly, 9) and yu(?/, 9) are the surface densities of these potentials. By the 
construction the heat potentials V{x^t) and W{x,t) vanish at t = 0. 

For bounded density ui^x^t) the heat potential of a single layer V{x,t) is continuous 
everywhere in M^, also on passing through the surface 5, and satisfies the homogeneous 
heat-conduction equation ()2.1()|1 outside S", i.e., it is parabolic outside S. The normal deriva- 
tives of V{x,t) have jump discontinuities on 5*. For continuous in 5* density v{x,t) these 
discontinuities are given by 

'dVix,t)\ dVix.t) 1 , , 



dn^ /. dn^ 2 

dV{x,t)\ dV{x,t) 1 , , ^ 

^ ' ' ^ ^ ' ' +-iy(x,t), xe S. (2.19) 



dn^ /e dn^ 2 

For bounded density /i(x, t) the heat potential of a double layer ()2.18|) is continuous 
everywhere outside of 5* (in Mf^yS") and in S. Outside S the potential W{x,t) is parabolic. 
On passing through S it has discontinuities. For continuous in S density /i(x, t) these 
discontinuities are given by {d = 3) 

Wi{x,t) = W{x,t) + ^fi{x,t), 

We{x,t) = W{x,t) -^fi{x,t), xeS. (2.20) 

The normal derivatives of the double-layer potential are continuous on passing through S. 

The employment of the Newtonian and heat potentials for transforming the boundary- 
value problems for Laplace and heat equations into the integral ones is based on the discon- 
tinuity properties on the boundary of the double layer potential and the normal derivatives 
of the single layer potential. Let us consider a simple example, namely, construction of the 
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solution to the Dirichlet problem for the heat equation in a compact domain D bounded by 
a smooth surface S: 

d \ 

— - a'^Aj u{x,t) = 0, xeD t>0, (2.21) 
u{x, 0) = 0, X G -D , 
u{x,t) = ip{x,t), X & S, t>0, 

where the function tp{x,t) specifies the temperature on the boundary 5* at different time 
instants t. We shall look for the solution u{x,t) in terms of the heat potential of a double 
layer (jTT^ 

u{x,t) = W{x,t) . (2.22) 
With account of Eq. ()2.20p we have on the boundary 5* 

\n{x,t) = -a^! dt' [ dS^'^^{x,t;x\t') fi{x',t') +^{x,t), x,x' e S . (2.23) 

2 Jo Js (y^x' 

Thus, the problem under consideration is reduced to the solution of the linear integral 
equation of the second kind. With respect to the variable t these equations are of the 
Volterra type and with respect to the spatial variables x and x' they are of the Fredholm 
type, the variables x and x' ranging on the boundary S. 

In an analogous way the integral equations for the Green's function of the heat equation 
can be deduced. Let us consider this technique in the case of the first boundary-value 
problem (the Dirichlet problem) for this equation. The Green function^ K{x, t; x', t') is 
specified by the following conditions;^ it should satisfy the homogeneous heat equation with 
respect to the first pair of its arguments 

^ a^A^. ) K{x, t; x' , t') = , (2.24) 



it should obey the inhomogeneous initial condition 

K{x,t;x',t) = 6{x,x'), t>0, (2.25) 

and the homogeneous boundary condition 

K{x,t;x',t') = 0, xeS. (2.26) 

We represent the Green function K{x,t; x' ,t') as the sum of a free propagator and the 
heat potential of a double layer 

K{x, t; x', t') = Ko{x, t; x' , t') + [ dO [ dSy ^{x, t; y, 6) /i(?/, 6; x\ t') . (2.27) 

Jt' Js (J^y 
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The right-hand side of this equation obviously satisfies Eq. ()2.24p . When t = t' the double 
layer potential in ()2.27p (the second term) vanishes. The free propagator in this formula Kq 
enables one to obey the initial condition (|2.25|) . The density of the double layer potential 
fi{y,6] x' ,t') is determined from the boundary condition (j2.26p . On substituting Eq. (j2.27p 
into (j2.26|) the following integral equation is obtained for the potential density /i 

- ^{x,t; x' ,t') = —KQ{x,t; x' ,t') — [ dO fdSy^--^{x,t;y,6)^{y,6;x\t'), xGS. 
2 Jt' Js '^^y 

(2.28) 

For Eq. ()2.28j) and consequently for Eq. ()2.27|1 the method of successive approximations 
can be developed, as Eq. ()2.28|1 is an integral equation of the second kind. It has been 
proved, that the series arising here is uniformly convergent.— The first terms of this series 
for Eq. ^(m\ are 

K{x,t;x',t') = Ko{x,t;x',t') + (-20^)^ / d9 dSy —^{x,t;y,9) Koiy,e;x',t') 

Jt' Js c'riy 

+ {-2a^f / d9 dSy—^{x,t-y,e) / d9, / dSy,—^{y,9-y,,ei)Ko{yi,e,;x',t') + .... 
Jt' Js '^'"■y Jt' J (^'^yi 

(2.29) 

Obviously, this series is a result of successive approximations applied to the integral equation 
for the complete propagator 

/"* r dK 

K{x,t;x',t') = Ko{x,t;x',t') -2a'^ d9 dSy —^{x,t;y,9) K{y,9;x',t') . (2.30) 

Jt' Js ' ^^y 

The series fl2.29p is nothing else as the multiple scattering expansion for the heat kernel 
in the problems under consideration . Thus we have derived this expansion in a rigorous 
way. 

By making use of the Laplace transform in Eq. (j2.29|l one can remove the integrations 
over the intermediate time variables ^'s 

K{x,x';p) = Ko{rxx';p) + {-'2a'^y / dSy Ki{rxy;p) Ko{ryx';p) 

Js 

+ {-2a^f / dSyKi{r^y;p) / dSy^ Ki{ryy^;p) Ko{ry^^r,p) + . . . , (2.31) 
Js Js 

where 

POD 

K{x,x';p) = / e-P^K{x,t;x\0)dt, (2.32) 
Jo 
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Koir,y;p) = / e-^'Koix,t;y,0)dt=—^Ko(^Vp)^ (2-33) 
K,{T^y-p)= re-^'|^(^>^;?/>0)rft=-^VP^i(— v^)> y^S. (2.34) 

^0 Tfd Qj 

In Eq. ()2.34p ip is the angle between the vectors r^y and rij^. The Laplace transforms Kq 
and Ki are calculated for d = 2. They are expressed in terms of the modified Bessel 
functions^^ Kq{z) and Ki{z). We hope that our notations will not lead to confusion because 
the free propagator Kq{x, t; x' , t') and the Bessel function Kq{z) have different number of the 
arguments. The series (j2.3H) is the perturbative solution to the following integral equation 

K{x,x';p) = Ko{r^^r,p) -2a'^ dSy Ki{r^y;p) K{ry^r,p) . (2.35) 

Js 

The series (j2.29j) or (j2.31|) contains complete information about the Green function (heat 
kernel) in the problem at hand. However, extracting it from here is not a simple task. 

The second term at the right hand side of Eq. (j2.31|) is responsible for one 'reflection' 
from the boundary (the Born approximation). Its contribution into the heat kernel can be 
expressed in terms of the confluent hypergeometric functioni^ Wajs- By making use of the 
convolution theorem for the Laplace transforn>2ii2i22iSi we obtain 

K('\x,x';p) = l^dSyCos^Ko (^v^) K, (^^) . (2.36) 

The inverse Laplace transform givea^^ 

- ^J^, exp {-^)W.. , ,2.37) 

where (f is the angle between the vector r^y and the inward directed normal to the boundary 
5* at the point y. 

The perturbation series ()2.31|) can be employed, for example, to find the asymptotic 
expansion of the heat kernel trace. 

III. COMPOUND MEDIA 

An important advantage of the heat potential technique for constructing the integral 
equations is the possibility of applying it to compound media. We show this by considering 
first the solution of the heat-conduction equation instead of the relevant Green's function. 
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Thus, in both the regions and D_ the heat equations are defined 
d \ 

— - a^^Aj u+{x,t) = Ttx{a+)u+{x,t) =0, x e D+ , (3.1) 
d \ 

a^_A] u^{x,t) = Tt^{a^)u4x,t) =0, x e (3.2) 



with the matching conditions at the interface S, namely, when crossing S the following 
quantities should be continuous: temperature 

u^{x,t) = u_{x,t), xeS (3.3) 

and heat current 

A,^±(£lf + A_^z(^ = 0. .eS, (3.4) 

on^[x) on_[x) 

where ?7,+ (x) and n^{x) are inward normals to the surface S at the point x for the regions 
Dj^ and -D_, respectively. These matching conditions imply, in particular, that there are no 
heat sources on S. The parameters a+, a_, A+, and A_ specify the material characteristics 
of the media. 

We shall look for the solution to this problem in terms of the heat potentials of single 
layer and double layer. Here the following feature proves to be important. If the solution 
u+{x,t) in the internal region Dj^ is represented as the heat potential of a single layer 

u+{x,t) = a\ f dt' [dS^'K^o^\x,t; x' ,t') u{x' ,t') , (3.5) 
Jo Js 

then the solution U- {x, t) in the external region D_ should be looked for in terms of the 
heat potential of a double layer 

U-(x,t) = a^_ / dt' dSx'T: — 7 — r (x, t; x', t') /x(x', t') , (3.6) 
Jo Js dn_{x') 

where K^'^^ and ^ are the fundamental solutions of the heat equations ()3.H) and ()3.2|) . 
which are defined by the formula ()2.12|) with a = a+ and a = a_, respectively. 

Substituting Eqs. ()3.5p and ()3.6|) into the first matching condition ()3.3p we obtain 

ft 



f dt' [ dS^^K^;^^ {x, t; x', t') z/(x', t') 
Jo Js 



^{x,t) + a^_ I dt' / dSx>- — '^—{x,t]x',t')^ix',t'), x, x' e S . (3.7) 



dn 
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The second matching condition ()3.4p results in another integral equation 

/■* r dK^'^^ 1 
A+a^ / dt' / dS,> ° (x, t- x', t') v{x', t') - - X+iy{x, t) 
Jo Js on^{x) Z 

+ \_al f dt' [ dS.,> \ , (a;, t- x'.t') t') = 0, x,x'eS. (3.8) 

Jq Js dn^{x)dn^{x') 

Thus the problem under consideration is reduced to the solution of the system of two 
linear integral equations of the second kind ()3.7|1 and ()3.8|1 . It is worth noting that we have 
obtained homogeneous equations, because there are no any heat sources in the problem 
under study. Hence, we are dealing here with the eigenfunctions only. 

Let us proceed to the Green's function K{x, t; x', t') in this problem. In what follows, it 
is convenient to represent this function in terms of the following four components depending 
on the range of the arguments: 

' K^^{x,t; x' ,t'), X, x' G -D+ , 
K+_{x,t]x',t'), xeD+,x'eD^, 



K{x,t;x',t') 



(3.9) 

K^+{x,t;x',t'), xeD^,x'eD+, 



{x,t;x',t'), x,x'eD_. 

The conditions, which specify the Green function K{x, t; x', t'), can be found in the follow- 
ing way. This function should provide the solution u{x,t) of the inhomogeneous boundary- 
value problem ()3.H) . ()3.2|) . ()3.3|) . and ()3.4|) with the heat source f{x,t) in the form 

u{x,t) = -e{t-t') [ dr [ K{x,t;^,T)f{^,T)d^, (3.10) 

Jt' JRa 

where 6{t — t') is the step function. The function u{x,t) will satisfy the inhomogeneous 
heat-conduction equation 

{Ttxia+) uix,t), X E D:, 
ftx{a^)u{x,t), xeD_, 

if the Green function K{x,t; x' ,t') obeys the corresponding homogeneous heat equations 
with respect to the first pair of its arguments 

Ttx{a+) K{x,t;x',t') = 0, x e D+ , 

ftx{a^)K{x,t;x',t') = 0, x e D_ (3.12) 
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and the inhomogeneous initial condition ()2.25p with respect to the both pairs of its argu- 
ments. In terms of the components ()3.9|) the initial condition ()2.25|) acquires the form 

K^+{x,t; x' ,t) = 6{x,x'), x, x' G -D+ , (3.13) 

K {x,t;x',t) = 6{x,x'), x,x'eD^, (3.14) 

K+_{x,t;x\t) = K_+{x,t;x\t) = 0. (3.15) 

On the right-hand side of the initial conditions ()3.15|) the delta function S{x, x') with 
x^x' & S is absent. Thus in our consideration we eliminate the treatment of heat sources at 
the interface. The point is such sources alter the matching conditions instead of the initial 
conditions ()3.4j) . In the next section the solution to the heat conduction equation, defined 
on a line, will be constructed for such a configuration by making use of the heat potential 
technique. 

The matching conditions ()3.3p and ()3.4|) are directly transformed into the conditions for 
the Green function K{x,t; x' ,t') with respect to the first pair of its arguments 

K+^{x,t; x' ,t') = K_^{x,t; x' ,t') , (3.16) 
' 'x,t;x',t') + A_ ^ \ {x,t;x',t') =0, (3.17) 



dn^{x) ' ' ' dn^{x) 

K+_{x,t]x',t') = K^-{x,t;x',t') , (3.18) 

dK dK 

^~ 'x,t;x',t') + A_- — ^{x,t;x',t') =0, (3.19) 



dn^ (x) 5n_ (x) 

X E S . 

It turns out that the heat equations ()3.12|) . the initial conditions ()3.13|) - ()3.15|) . and the 
matching conditions ()3.16|) - ()3.19|) are enough for construction of the Green function in the 
case of compound media in a unique way. 

We shall look for the components of the Green function (j3.9|) in terms of the heat po- 
tentials of a single and double layers with respect to the first pair of their arguments, the 
components and being expressed through the heat potentials of single layers and 

the components K and K^^ through the heat potentials of double layers. In order to take 

into account the inhomogeneous initial conditions ()3.13p and ()3.14j) for the components 

and K we add to the chosen heat potentials (nonsingular part of the Green function) the 

free propagator K^'^^ or Kq~'^ (singular part of this function) 

K++{x,t;x',t') = K^'^\x,t;x',t') + al [ dO [dSy K^'^\x,t;y,e) u++{y,e;x',t'),{3.20) 

Jt' Js 
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K+_{x,t;x',t') = al [ de [dSy Ki^'^\x,t;y,e) u+^{y,e; x' ,t') , (3.21) 

Jt' Js 

ft f f)K^~^ 

K(x,t;x',t') = K^o~\x,t;x',t') + al dO dSy ^ ° A x,t;y,e) u^^{y, 9; x\t' )(,?,. 22) 

Jt' Js ■ dn^{y} 

K_+{x,t-x\t') = al j^dO jdSy^-^^{x,t-y,e)u.+ {y,e-x\t'). (3.23) 

The matching conditions ()3.16p and ()3.17|) give 

K^^\x,t;x,t')+al [ dO [dSy K^'^\x,t;y,e) u++{y,e; x' ,t') 
Jt' Js 

1 /"* /" dK^^^ 

-z/_+(x,t;x',t') + a'- / de j dSy -^^^ix,t;y,e) iy_+{y,9;x\t') , (3.24) 



y 

t' Js 



dK ' f f dK^ ' 

X+- — ^{x,t;x',t') + X+al / de dSy - — ^(x, t; y, 6*) z/++(?/, 6*; x', t') 
dn+{x) ^ Jt, Js ' dn+{x) 



-^K'^++iy,0;x',t') + \_al J de jdSy ^^ ^^^^^ ^^^ {x,t;y,e) iy_+{y,e;x',t') = , 

xeS. (3.25) 
In the same way we deduce from ()3.18|) and ()3.19p 

/ de [dSyKl+\x,t;y,e)u+_iy,e;x',t') 
Jt' Js 



1 f dK^^^ 

= K^o'\x, t; x', t') + -u—{x, t- x\ f) + al j^ de JdSy g^^i^^ ^) '"—(y, x\ t') , 

(3.26) 

1 f dK^~^^ 

-- K'^+-{x,t;x',t') + X+al_ / de dSy - — ^{x,t;x' ,t') u+^{y,e;x',t') 

2 ^ Jt' Js dn+{x) 

ft f d'^K^^^ dK^~^ 

+X_al de dSy- — — ^ — —{x,t;y,e)u__{y,e;x',t') + X_-—^{x,t;x',t') = 0, 



I, Js dn^(x)dn^(y) dn^{x) 

xeS. (3.27) 

The sets of integral equations of the second kind ()3.24p . ()3.25p and ()3.26p . ()3.27|) define 
the heat kernel for compound media in full. With respect to spatial variables these equations 
are of Fredholm type while regarding time variable they are of Volterra type. It is essential 
that the integration over the spatial variables is restricted by the interface S only. Hence the 
dimension of the initial problem is reduced by 1. By making use of the Laplace transform 
one can remove the integration over the time variables in Eqs. ()3.24|) - ()3.27p as it has been 
done in Sec. II. 
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Obviously the integral equations for the heat kernel derived here can be also applied when 
the surface S divides the same medium into the regions and D^, i.e. when the constants 
and a?_ equal. 

For constructing the solutions to the integral equations (j3.24j) - (j3.27|) the perturbation 
theory can be employed (see Sec. II). The expansion parameters in this case prove to be the 
constants and A±. The perturbation series generated here are nothing else as the multiple 
scattering expansion for the heat kernel . Thus we have proposed a rigorous derivation 
of these expansions both for homogeneous media and compact regions and for compound 
media. 



IV. ASYMPTOTIC EXPANSION OF HEAT KERNEL FROM PERTURBATION 
SERIES 

In practical applications, especially in QFT, the asymptotic expansion of the integrated 
heat kernel when t — > +0 proves to be important . ^i^i^i^ It has the form 

/oo 
dxK{x,t;x,0) = (Ant)-'^/'^ ^ T/^E^/a + ES. (4.1) 
n=0,l,2,... 

In this expansion d is the dimension of the configuration space and ES stands for the ex- 
ponentially small corrections as t — > +0. We show how to derive this expansion proceeding 
from the perturbation series for the integral equations (|3.24|) - (|3.27|) . The functions K__^_ 
and do not contribute to integrated heat kernel (see subsection IV C|) , thus we have to 

consider only 7^++ and K 

For our purposes it is convenient to use such coordinates that in the vicinity of the surface 
5* the metric is Qij dx^ dx^ = (dx^)"^ + Qat dx'^ dx^ where x^ is a coordinate on the normal to 
5, x^=0 on S. In view of the exact form of the free propagator ()2.12j) . one can infer that 
in each term of the perturbation series for integral equations fl3.24|l - ()3.27j) the power in 
t contributions are given only at the following conditions: when evaluating the heat kernel 
trace, the integration over dx should be spread over the region immediately adjacent to the 
boundary S and in the course of the multiple integration over the boundary 5* the respective 
distances Tyyi should be also small. Therefore in the vicinity of S we may replace the squared 
distance {x — z)"^ by several terms of its expansion in powers of the corresponding geodesic 
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distance a on the surface 5" 



[X — z] 



(xs - + cr^l-{x3 + y3)ki + X3Z3{kl + kl)} 



-{2Z3 + X3)k[ + X3Z3{kik[ + k2k2, 



+ ... 



(4.2) 



k[ 



dki 
da' 



dk2 
da 



The surface area element is dSz = (l — ^^-Rab^"^^ cr^ + . . .) adadQ, Q parameterizes a 
unit sphere, Lab is the second fundamental form on S, Rab is intrinsic Ricci curvature, ^ is 
a unit tangent vector at x to the geodesies with the length a joining 2; to x on (see, for 
example, Ref. |23|). 

Here we present the first three terms of the perturbation series under consideration when 
(rf = 3) 



t-i s 



(47ra2.)3/2 



(47ra^)3/2 



Stt A+ a\ 



f 







vr3/2a+ y/^^^ + 28vr75 



1 A_rV2 



32 vr 



87r3/2A+ a+ js 
A_ (a+ + 2a_) 



d5 



5 (LlY + L'^Ll - -R: 



A+ (a+ + a_ 



5 



-(L:)^ + 4L^L»--i?," 



+ 



1 A_ 



1 



A 1 (a_L + a_ 



9a^ 



3a^ 



3 + 4a, 4al 



35n llo r\2 

— a H a a+ + 2a ; a_ 

12 " 4 " + + 

^rfS [(L^)2 + 2LX]|+... 



To obtain i^i^l(t) and iri^l(t) one should replace a+ ^ a_, A+ ^ A_ 
asymptotics of the subsequent terms of perturbation series may be found in a similar way. 
After that all factors appearing with the same powers of t are added up to give the heat 
kernel coefficients. The latter are expressed through the integrals of the surface geometric in- 
variants. The asymptotics ()4.3j) were presented in Ref. without considering the derivation 
of the relevant integral equations. 



(4.3) 
-Ll The 



V. HEAT KERNEL ON A LINE 



In this section we demonstrate the efficiency of our approach based on integral equations 
for constructing the heat kernel on a line. In this case the interface between the media 
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reduces to a point. As a result we are dealing with the Volterra integral equations in respect 
of one (time) variable. These equations are of a special type (Abel equations), and their 
solutions can be found in an exact form. 



A. Homogeneous media with gluing conditions 

By making use of the heat potential technique we construct here, in an exact form, the 
heat kernel K{x, y; t) for the Laplace operator on an infinite line for homogeneous medium 
with a nonstandard gluing conditions at the origin (these conditions will be specified below) . 
From the physical point of view K{x, y; t) is the temperature at the point x which is generated 
by a unit instantaneous heat source placed at the point y at the moment t = 0. 

As in previous sections we first formulate the conditions that define the heat kernel in the 
problem under consideration. With respect to the first argument K{x,y]t) should satisfy 
the one-dimensional heat-conduction equation 

d 92 \ 

— - — 1 K{x, y;t) = 0, t > 0, -00 < x < 00 (5.1) 
and special conditions at the interface x = 



l-^^K{x,y;t) 



lK{-0,y;t) = l-'K{+0,y;t), (5.2) 

(5.3) 



d 

= l—K{x,y;t) 
x=-o 9x 



x=+0 

where / is a dimensionless parameter. We shall refer to these conditions as to gluing ones. 
In the Appendix A it is shown that such conditions lead to a selfadjoint spectral problem 
for the Laplace operator in any dimension. The initial condition for K{x, y; t) involves its 
both space arguments 

K{x,y;0) = 6{x,y). (5.4) 

We shall seek for K{x, y; t) in terms of free heat kernel and single layer heat potentials.— 
The solution is decomposed in four components related to different positions of the heat 
source and the observer 

K^+{x,y;t) = / dr Ko{x,0;t — t) ai{T,y), x < 0, y > , (5.5) 
Jo 

K++{x,y;t) = Ko{x,y;t) + dr Ko{x,0;t - t) a2{T,y), x,y>0, (5.6) 

Jo 



17 



K+_{x,y;t) = / dr Ko{x,0;t - t) a3{T,y), x > 0, y < , (5.7) 

^0 

K__{x,y;t) = Ko{x,y;t) + dr Ko{x,0;t - t) a4{T,y), x,y<0, (5.8) 

Jo 



where Ko{x,y;t) is the free heat kernel (propagator) on an infinite fine 

1 (^-yf 



Ko{x,y]t) 



e 4t 



(5.9) 



First we substitute Eqs. ()5.5|1 and ()5.(i|l into the gluing conditions ()5.2|1 and ()5.H|1 . Then 
we take into account that the single layer potential changes smoothly across the boundary, 
while the normal derivative of this potential undergoes a jump 

1 



K_4-0,y;t) = 



2v^ Jo \/t-r 



K++i+0,y;t) = KoiO,y;t) + -^ f dr^^^ 

•^V^ Jo Vt-T 



d_ 
dx 
d_ 
dx 



K.+ {x,y;t) 
K++{x,y;t) 



x=-0 



x=+0 



ai{y,t) 



d 



2 "2(y,t) + ■^Ko{x,y;t) 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



x=0 



Inserting Eqs. (j5.1Up and (j5.11|) into Eq. (j5.2p one obtains the Abel integral equation 



dr 



Jo \/t -T 



[/^ai(y,r) - a2{y]T)] = Ko{0,y;t) 



with the solution^ 



ai{y,t)-l ^«2(z/;t) = TTT^ / -^0(0^^) 



P^/tt dt Jo \Jt — T 
The substitution of Eqs. (j5.12|) and (j5.13p into Eq. (j5.3j) gives 

ai{y,t) = -f a2{y,t) + ^-^^e 4i 

From (fH?TK|) and (fCT)|) it follows that 



1 y _r_ 



ai{y,t) 



/4 + 1 0Ft3/2 



And finally 



K_+{x,y;t) 
K++{x,y;t) 



e 4t, a2{y,t) 



e 4t 



/4-1 



/4 + 1 2v^t3/2 



e 4t, 



1 _(£zl)! Z^-1 1 

e 4t + — e 4i 



2V7rt 



+ 1 2x/H 



(5.14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 
(5.19) 
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1 _(x-yf 

e 4t 



In a similar way one gets 

K+-{x,y;t) 
K__{x,y;t) 

These formulas are is in a complete agreement with the resu 



1 i^-y)^ 



/^-l 1 



it — 



+ 1 2v^ 



4t 



(5.20) 
(5.21) 



of the Lemma 5.2 argued in Ref. 



25 and Lemma 4.1 from Ref. 



t of a combined employment 
29 



K++{x,y]t) = cos"^ 9KN{x,y;t) + sin^ 9KD{x,y;t) , 

K {x,y]t) = siia^ 9Kiy{x,y;t) + cos^ 9KD{x,y;t) , 

K+_{x,y]t) = sm9 cos9[KN{x,y;t) - KD{x,y;t)] , 



(5.22) 
(5.23) 
(5.24) 



where 



cos^^ 



I' 



sm^9 



1 



i^ + r i^ + i' 

and K£){x,y;t) and K]sf{x,y;t) are the heat kernels for Dirichlet and Neumann boundary 
conditions, respectively. 



B. Dielectric-like conditions on a line 

We construct here the heat kernel for the Laplace operator defined on an infinite line 
with dielectric-like matching conditions at the origin a; = 0. The heat kernel is defined by 
the heat conduction equation 



d 



dt 



dx'^ 



— a {x)Tr^ K{x^ y;t) = 0, t > 0, — oo < x < +oo, x . 



(5.25) 



where 



a^x) 



X < 0, 
X > 0, 



and being positive constants. At the interface of dielectric media the matching 
conditions 



K{-0,y;t) = K{+0,y;t) 



d 



X-^K{x,y;t) 



x=-0 



d 



(5.26) 
(5.27) 



x=+0 



should be met. As usual, the initial condition for K{x,y]t) is given by ()5.4|) . 
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We call the boundary conditions ()5.26p and ()5.27|) the dielectric-like conditions. The use 
of this term requires some explanations. When two dielectric media and possessing 
different characteristics are separated by the interface S of an arbitrary form then in Maxwell 
theory we have on the surface 5* the set of coupled boundary conditions involving all the 
components of the electromagnetic potential yl^()f:, x), /x = 0, . . . , c?. If one disregards the 
vector character of the electromagnetic field and confine oneself to oscillations described by 
a sole scalar potential (for example, sound waves which are described by a scalar velocity 
potential) then at the interface between different media the conditions ()5.26|) and ()5.27|) 
should be satisfied. In other words, these boundary conditions hold in the theory of scalar 
'photons' in compound media. 

Again we seek for the solution in terms of a relevant free propagator and single layer heat 
potentials: 



i\:_+(x, y] t) = ai I drKoix, 0; ai(t - r)) /?i(r, y), x <0,y>0, (5.28) 
Jo 

K++{x, y; t) = Ko{x, y; a^t) + a\ j drKoix, 0; a+(t - r)) (32{r, y), x,y > 0, (5.29) 

Jo 

K+_{x, y- t) = al [ drKoix, 0; a^(t - r)) P^ir, y), x>0,y<0, (5.30) 
Jo 

K {x,y]t) = Ko{x,y]a^_t) + o?_ I drKo{x,Q]a?_{t — r))j3/^{r,y), x,y<0. (5.31) 



First we insert Eqs. fl5.28|) and ()5.29|) into the matching conditions ()5.26p and ()5.27p . 
The single layer potential changes smoothly across the boundary, while its normal derivative 
undergoes a jump 

i^-.(-0,.;t) = ^r^-^> (5-32) 



2V^Jo 



K^4+0,y;t) = Ko{0,y;alt) + ^ [ dr^^^ , (5.33) 

2vvr Jo y/t-r 



d , , 
—K_+{x,y;t) 

-^K++{x,y;t) 



= -lf3i{y,t), (5.34) 
= -\i52{y,t)^ ^Ko{x,y;a\i) . (5.35) 



x=0 

Substituting Eqs. ()5.32|) and ()5.33|) into Eq. ()5.26|) one obtains the Abel integral equation 

\ dr , ^ , . „ 2_ 



Jo \Jt- r 



[a_/3i(y, r) - a+/?2(l/; t)\ = Ko{0, y; a,t) (5.36) 
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with the solution 

ft 



2, d 1 diT 

a-pi{y,t) - a+p2{y;t) = Ko{0,y;alT) 



^ ^ expf--4T). (5.37) 



The substitution of Eqs. ()5.34|) and ()5.35|) into Eq. ()5.27|) gives 

A_A(y,t) = -X^p,(y,t) + X^^^^e'^' . (5.38) 
From Eqs. (jOTjl and it follows that 



pjy^t) = °+ , ^ ,/o e /?2(l/,t) = 4± ^ e (5.39) 



y 



X+ X- 



And finally 

A+ _(£-y_a_/a+)f_ 

^-+(^,2/;^) = ^'^^-^ , (5.40) 



H A/Trai t 

a+ a_ V + 



In a similar way one gets 



2 (x-ya+/a-)^ 



^+-(^,2/;t) = ^T^^^e , (5.42) 
— - \ rn a t 

^"(^'^'^) = 7^^^"'"-* +^r-^7^^^"^- (5-43) 

2\/Tcatt h— 2A/7ralt 

V a_ a+ V 

The solution obtained here exactly reproduces the results obtained in this problem by other 
methods.i2»2I 



C. (5-Like heat source at the interface 



In preceding considerations we excluded the configuration when the 5-like heat source is 
placed at the interface of the media y = 0. For completeness we have to check whether this 
configuration contributes to the trace K{x, x; t) dx. To this end we use the approach of 
Ref. |2S!. The idea is to modify the boundary conditions ()5.26|1 and ()5.27|) so that they allow 
for the heat source placed at the interface 

ir_(-0;t)=ir+(+0;t), (5.44) 
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x=-0 



d 



QS{t) 



(5.45) 



x=+0 



The condition (j5.45p means that the heat Q instantly generated by the source is divided 
into two flows which are proportional to A_ and A+. We represent the solution in the form 

ft 



K_{x;t) = at / dTKo{x,0]ai{t - t)) f3i{T), x<0, 



K+{x]t) = a\, I dTKo{x,0;a^{t - t)) (32{t), x>0. 



(5.46) 
(5.47) 



The functions /3i(r) and /52(t) can be determined by making use of the Laplace transform. 
We denote by K{x; s) the transform of K{x; t), i.e. 

POO 

K{x;p) = / dte^'P' K{x]t). 



Transforming ()5.4fi|l and ()5.47|1 we obtain 



K_{x;p) = a_/5_(p) V'vtTp exp(-x^/a„), 



K+{x;p) = a+p+{p)^/7^ exp{-xy/p/a+). 
Then Eq. (I5.44|) leads to the relation between the Laplace transforms /5_ and 

a-P-{p) - a+f3+{p) = 0. 
The substitution of Eqs. ()5.46|) and ()5.47|) into Eq. ()5.45|) gives 



A_/5_(t) + A+/5+(t) = 2Q5(t). 



After the Laplace transform one arrives at the second relation between (3i and (32 



X,p.{p)+X2P+{p) = 2Q6{p). 



(5.48) 
(5.49) 



(5.50) 



(5.51) 



(5.52) 



How to apply the Laplace transform to the singular 6 function can be found in appropriate 
handbooks . The essence of the matter comes to defining the integration rule 



Jo 



f{t)6{t)dt = f{0), 



whence it follows in particular 



S{p) = / e-P'6{t) dt=l. 
Jo 
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In the problems treated by the integral Laplace transform the semiaxis t > (or t > to) is 
usually considered. Therefore the 5-function should be defined here in a nonsymmetric way, 
for example, as the limit when e — > +0 of the function 

0, t <0, t > 

<t < e. 

The solution of the system ()5.5()j] . ()5.52|1 is 



2a,Qm 2a_QSip) 

A+ a_ + A- a+ A+ a_ + A_ a+ 



By the inverse Laplace transform we find from ()5.53|1 



A+a_ + A_a+ A+ a_ + A_ a_|_ 



Having inserted ()5.54j] into ()5.4(i|l and ()5.47|1 we derive 



K4x;t) = -^e ^^'^ ^ , x<0, (5.55) 



a+ 



1 - O 

K+{x]t) = -y= e ^ ^ ^ ^ , x>0. (5.56) 



TTt ^ + 

a+ a_ 

These formulae show that for t > the heat kernel is finite notwithstanding the S-\ike heat 
source situated at the interface between two media. Therefore the neighborhood of the point 
a; = gives no contribution to the heat kernel trace K{x, x; t) dx. 



VI. CONCLUSION 



For a broad set of boundary conditions the finding of the heat kernel is reduced to the 
solution of integral equations defined on the boundary (or at the interface) of the manifolds. 
As a result the dimension of the initial problem is brought down by 1. Remarkably this 
technique is applicable to compound media where the standard methods for the investigation 
of heat kernel do not work because in this case the principal part of the elliptic operator in 
question is not smooth. 

The perturbation series for the integral equations derived are nothing else as the multiple 
scattering expansions for the relevant heat kernels. Thus a rigorous derivation of these 
expansions both for homogeneous media and compact regions and for compound media has 
been done. 
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The efficiency of this approach is convincingly demonstrated by constructing, in an exact 
form, the heat kernel on an infinite line with diverse matching conditions and by deriving 
the first terms of the asymptotic expansion for integrated heat kernel in the case of three 
dimensional compound media. 
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APPENDIX A: SELF-ADJOINTNESS OF BOUNDARY VALUE PROBLEMS 
FOR COMPOUND MEDIA 

Often it is helpful to use the spectral representation for the Green's function of the heat 
conduction equation ()2.10p 

K{x,t;x',0) = Y,e^^''h{x)mx'), (Al) 

k 

where fk{x) are the eigenfunctions of the spectral problem at hand 

- a^Afkix) = LUkfkix), xeD (A2) 

obeying the relevant boundary conditions on 5*. Apparently, the representation (jAljl is well 
defined if the spectral problem ()A2|) is hermitian and positive definite. In this connection it is 
worth elucidating the boundary conditions for compound media that lead to the self-adjoint 
spectral problem (see Eqs. p.ip - ()3.4p ). In this case we have instead of ()A2|) 



- a^{x) Afk{x) = ujkfk{x) , (A3) 

where 



a 



+ ; X G -D+ , 



a\x) = i I (A4) 
ai, x G , 

a\ and being constants. At the interface S the natural modes fk{x) obey the dielectric 
like conditions 

fk+{.x) = fk-{x) , 
24 



^ dfk+{x) ^ ^ dfk-{x) _ ^ 
9n+(x) dn^{x) ' 

xeS. (A5) 

Here we are using the same notations as in Eqs. ()3.1|) - ()3.4|) . 

Since the outer region D_ is not bounded the spectrum uj^ is continuous. For example, we 
can assume that at the spatial infinity |x| — oo the natural modes fk{x) satisfy the scattering 
problem conditions and decrease sufficiently fast. The explicit form of these conditions will 
not be needed below. 

The differential operator in Eq. ()A3jl apparently coincides with its adjoint 

[-a\x)A]]=\ ^ ^ (A6) 

y -a_A, X e D_. 

The boundary value problem with the operator — a^(x)A will be self-adjoint when the 
integral for two sufficiently smooth functions u{x) and v{x) 

dxa'^{x){vAu-uAv) (A7) 

D+UD- 

vanishes. Applying the Green integral formula for the domains and D_ separately we 
obtain 



dS 



s 



2 / dv+ du+\ 2 / dv^ du- 
I u+— Vj^— — I + a_ I u^- f_ 



(A8) 



The functions u{x) and v{x) are assumed to diminish at the infinity in such a way that the 
region |x| oo does not contribute to the integral ()A8|1 . The integral / is equal to zero, for 
example, for the following conditions at the interface S 

Uj^{x) = u^{x) , 

r, du+(x) r, du^ix) ^ , 

al-^ + al—^ = {), xeS. A9 
on+{xj on^[x) 

The function v{x) should satisfy the same boundary conditions on S. Thus the dielectric-like 
conditions ()A5|1 lead to self-adjoint boundary value problem if 

A+ _ A- 

Of course, conditions ()A9|) do not exhaust all the cases when the boundary value problem 
under consideration is self-adjoint. Let the surface S divides the same medium into the 
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domains and D-, i.e. a\ = a^_. We get self-adjoint spectral problem if impose at the 
interface 5* the following gluing conditions 

/ u+{x) = l~^u^{x) , 

= -^S, (AlO) 

where I is a dimensionless constant. The one-dimensional version of this problem has been 
considered in Sec. IVl 
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